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Abstract 

We investigate one-flavor QCD with an additional chiral scalar field. For a large 
domain in the space of coupling constants, this model belongs to the same universality 
class as QCD, and the effects of the scalar become unobservable. This is connected 
to a "bound-state fixed point" of the renormalization flow for which all memory of 
the microscopic scalar interactions is lost. The QCD domain includes a microscopic 
scalar potential with minima at nonzero field. On the other hand, for a scalar mass 
term below a critical value m^, the universality class is characterized by per- 
turbative spontaneous chiral symmetry breaking which renders the quarks massive. 
Our renormalization group analysis shows how this universality class is continuously 
connected with the QCD universality class. 



1 Introduction 

Universality of QCD means that predictions are independent of the details of the micro- 
scopic interactions. This is crucial for predictivity, since the precise form of the funda- 
mental interactions at very short distance scales is not known. In a large parameter space 
characterizing possible fundamental interactions, the QCD universality class corresponds, 
however, only to a certain domain. For other domains in parameter space, the color sym- 
metry may be "spontaneously broken" by the Higgs mechanism, or all quarks may acquire 
a large mass due to spontaneous chiral symmetry breaking. We are interested here in 
the transition from one domain to another and in the question of what happens at the 
boundary of the "QCD domain" . 

Looking at QCD from a microscopic scale - say a unification scale lO^^GeV - its uni- 
versality class is characterized by eight massless gluons and a certain number of massless 
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fermions. Perturbatively, the masses are protected by the gauge symmetry and chiral sym- 
metries. At a much smaller scale around IGeV, nonperturbative effects induce masses for 
all physical particles. In particular, the fermions become massive owing to chiral symmetry 
breaking (xSB). This may be described by a nonzero expectation value a ~ (V'V') of ^ 
"composite" scalar field. In order to keep the discussion simple, we concentrate here on 
the case of one quark flavor - generalizations to several flavors are straightforward. 

Let us now consider a class of microscopic theories with a complex fundamental "chiral 
scalar field" (f) which has the same transformation properties as ipip and a classical potential 

^ = mV0+^A^(0»2. (1) 

The symmetries also allow for a Yukawa coupling between and the quarks. For nonzero 
(0), the chiral symmetry is broken and the quarks become massive. In the case of large 
enough positive im? (in units of some unification scale, say lO^^GeV), the scalar field is 
super-heavy and decouples from the low-energy theory. This range of im? obviously cor- 
responds to the universality class of QCD. All effects of the scalar field are suppressed by 
p^/m^, with p a characteristic momentum. For QCD predictions, they can be completely 
ignored. 

On the other hand, for large enough negative m?, we expect the perturbative picture of 
spontaneous symmetry breaking to hold. The scalar field gets a vacuum expectation value 
(VEV) 

(0) = a=|m^/A^,R|i/2^ (2) 

with mp{ and A(^^r related to m and by renormalization corrections. Both cr and the 
quark masses are of the order of the unification scale in this domain. The universality class 
now corresponds to gluodynamics without light quarks. In the chiral limit of a vanishing 
current quark mass, spontaneous xSB also generates a very light pseudo-Goldstone boson 
in addition to the gluonic degrees of freedom. 

Varying the microscopic scalar mass term m? from large negative to large positive values 
should lead us from the universality class with perturbative spontaneous chiral symmetry 
breaking (P^SB ) to the universality class of one-flavor QCD. One of the aims of this note 
is to understand the qualitative features of this transition in the vicinity of a critical value 
m^. This is clearly a nonperturbative problem, since on the QCD side of the transition the 
effective gauge coupling grows large. 

Our investigation is based on a nonperturbative flow equation which is obtained by a 
truncation of the exact renormalization group equation for the effective average action pQ . 
A crucial ingredient is the "bosonization" of effective multi-fermion interactions at every 
scale j2]. This provides for a description of fundamental scalar fields and bound states 
in a unified framework. A theoretical method with this feature is actually required for 
our problem, since the scalar quark- ant iquark bound states in the QCD description (e.g., 
the pseudo-Goldstone eta meson and the sigma meson) are expected to become associated 
with the fundamental scalar in the P^SB description. In this framework, we see also how 
one relevant and two marginal parameters in the PxSB universality class, namely the ones 
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corresponding to the mass and quartic self-interaction of the scalar field and the Yukawa 
coupling, become irrelevant for the QCD universality class. 

This remarkable change of the number of relevant parameters at the transition between 
the two universality classes is connected with the appearance of a bound-state fixed point 
for the flow of the scalar mass and self-interaction in the range of microscopic parameters 
corresponding to QCD. This bound-state fixed point is infrared attractive for all couplings 
except for the gauge coupling. Under the influence of this fixed point, all memory of the 
details of the microscopic interactions in the scalar sector is lost. This is exactly what is 
required for the QCD universality class which has the gauge coupling as the only marginal 
parameter (for a massless quark). In order to see the appearance of the bound state, it is 
crucial to re-incorporate the effective multi-fermion interactions generated by the flow into 
the effective bosonic interactions. This avoids an unwanted redundancy of the description. 
It also solves an old problem in the investigation of gauged Nambu-Jona-Lasinio models 
[3]; namely, how the presence of apparent relevant parameters in a too naive treatment of 
these models can be reconciled with QCD, where no such relevant parameters are present. 
In our approach, the flow towards the bound-state fixed point solves this generic problem. 

As a result of our investigation, we find a qualitatively convincing picture of the transi- 
tion between the two universality classes investigated. We have kept the truncation simple 
in order to illustrate the change in the number of relevant and marginal parameters in 
a simple way. The price to be paid is a limited accuracy in the quantitative description 
for parameter regions where the effective gauge coupling grows large. In our setting, this 
concerns primarily the quantitative details of the flow of the instanton-mediated interac- 
tions and the running of the strong gauge coupling. We emphasize that the qualitative 
picture does not require a detailed understanding of strong interactions in the momentum 
range where the gauge coupling is large. All decisive features are determined by the flow 
in a momentum range substantially larger than IGeV. In the same spirit, we also have 
neglected other effective bosonic degrees of freedom which may correspond to additional 
bound states. We keep only the composite scalars and the gluons. When we proceed with 
our analysis to the strongly coupled gauge sector, we do not attempt to compute the glu- 
odynamics, but simply model the strong interactions with an increasing gauge coupling; 
for the latter, we use various examples discussed in appendix^ We do not claim that our 
truncation of the gauge sector is sufficient in order to establish chiral symmetry breaking 
in QCD. A much more elaborate analysis would be needed for this purpose. We rather 
take the spontaneous symmetry breaking in the QCD universality class as a fact (estab- 
lished by other methods and observation). We only require that a reasonable truncation 
should describe chiral symmetry breaking. Beyond this, the details of the truncation in 
the gauge sector are not relevant for our discussion of universality classes. Despite these 
shortcomings, we expect that our quantitative results describe the right order of magnitude 
of one-flavor QCD. An impression of the size of uncertainties can be gained from Table. ^ 
in appendix lEl 

In order to illustrate our points, we compute the scalar condensate, i.e., the renormalized 
minimum of the effective potential, ctr = Z^l^oP, for a broad range of initial scalar 
mass values m^. We note that ctr, is directly connected with the decay constant of the 



3 




ml m\/k^ ml m\/k^ 



Figure 1: Renormalized scalar vacuum expectation value versus the initial condition 
for the scalar mass at the UV cutoff, fh\. Left panel: logarithmic plot for the parameters 
= 1, /i2 = 1, = 100 at A = lO^^GeV resulting in ml ^ -O.SSA^. Right panel: 
linear plot for fi^ = 0.1, = 1; for m\ < ml ~ — 0.0043A2, the linear dependence 
0"!^ ~ ~{m\ — ml) as expected from perturbation theory is confirmed. 



eta meson and sets the scale for the quark mass generated by xSB . We first neglect 
the anomalous Ua(1) violating contributions from instanton effects which only affect the 
physics at scales around IGeV. (They will be considered in Sect. El) We parametrize the 
microscopic interactions by the initial values of the renormalization flow at a GUT-like scale 
A = lO^^GeV. As can be read off from Fig. ^ a critical mass ml exists. For initial scalar 
masses below this critical mass, m\ < rhl, the naive expectation is fulfilled, and we find 
scalar condensates of the order of the cutoff, ctr, ~ 10^^ . . . lO^^GeV. It is remarkable that 
the value of the critical mass is negative and typically of the order of the cutoff or only a 
few orders of magnitude below the cutoff; for example, we find rhl ^ — O.SSA^ for the initial 
values = 1 and = 100 at A = 10^^ (Fig.^(left panel)). For a perturbatively accessible 
set of initial parameters h'^ = 0.1 and A^ = 1 at A = 10^^, we find rhl ^ — 0.0043A2 (Fig. [T] 
(right panel)). In the latter case, we find a linear dependence of the condensate on the 
mass parameter, cx^ ~ ~(^a ~ ^1): expected from perturbation theory (cf. Eq. Q). 

However, for initial scalar masses above this critical mass, m\ > ml, the scalar con- 
densate is 16 orders of magnitude smaller (not visible in the linear plot in Fig. ^ (right 
panel)). In this case, symmetry breaking is triggered by the fermion and gauge sectors and 
not by the scalar sector, i.e., ctr is roughly of the order of Aqcd- Therefore, even if we start 
the flow deep in the broken regime with m\ < but above the critical mass, the scalar 
fluctuations drive the system first into the symmetric regime where it will be attracted by 
the same IR fixed point as a QCD-like system. It should be stressed that no fine-tuning of 
the initial parameters is needed, neither to put the system into the domain of attraction of 
the QCD universality class nor to separate the UV scale from the scale of chiral symmetry 
breaking. 

Only for m\ < ml is the effective coupling between the scalars and the fermions strong 
enough to induce P^SB with a magnitude determined by the initial parameters of the 
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scalar sector. In this case, we would have to fine-tune the initial condition for Th\ to lie 
extremely closely to m^, if we wanted to separate the UV scale from the scale of chiral 
symmetry breaking. This is the famous naturalness problem which is generic for models 
involving a fundamental scalar. Of course, theories without fundamental scalars such as 
QCD do not have this problem, although effective scalar degrees of freedom such as bound 
states can occur at low energies. It is one of our main observations that the mechanism 
of how "QCD-like" theories circumvent the naturalness problem can also be applied to 
models with a fundamental scalar. 

The details of our study are organized as follows: in Sect. El we introduce the class of 
models containing one-flavor QCD and derive the flow equations for a qualitatively reliable 
truncation including "bosonization at all scales" . Section El is devoted to a discussion of 
the bound-state flxed point which governs the flow of the QCD domain for weak gauge 
coupling. In Sect.|31 we analyze the universal features of the QCD domain numerically and 
give estimates of IR observables in the nonperturbative strong-coupling regime. Instanton- 
mediated interactions are included in Sect. where we also describe the fate of the pseudo- 
Goldstone boson. 

2 Flow equations 

QCD with one massless Dirac fermion flavor coupled to an SU(A'^c) gauge fleld is charac- 
terized by the classical (or bare) action 

Sqcd = j d'xi, i]p[A] i, + (3) 

where = ^^J^S^^ — igT^^Al^^, and Ta denotes the (hermitean) generators of the gauge 

group in the fundamental representation. In this work, we embed one-flavor QCD in a 
larger class of chirally invariant theories including a color-singlet scalar fleld. For this, we 
consider the action 

r = y|z^^i^^ + ^[(^^)2-(^^5^)2] 

+ (4) 

which represents a simple truncation of the space of action functionals and serves as the 
basis of our approximations. Here we have used the shorthand (i/jip) = i^^i/Ji for the color 
indices. We included a background gauge flxing term with parameter ^, and = A^ + a^, 
An being the background and the fluctuation fleld, = D^lA]. Furthermore, we do 
not display the ghost sector for simplicity. Equation (jH) reduces to one-flavor QCD if 
we set the four-fermion and the Yukawa interaction equal to zero, = h = 0, let the 
scalar fleld be auxiliary, = 0, and set Zp = 1 = Z^ (the scalar potential is of no 
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importance then). Furthermore, there is a redundancy in Eq. Q: we can compensate for 
a shift in A^. by readjusting the Yukawa couphng and the scalar potential corresponding 
to a Hubbard-Stratonovich transformation (partial bosonization). But apart from this 
redundancy, which will be removed later on by "re-bosonization" , different initial values 
for the various parameters in Eq. generally correspond to different quantum theories. 
Some of these theories will belong to the same universality class sharing the same low- 
energy properties, which makes them indistinguishable from a low-energy physicist's point 
of view. 

We analyze this class of theories in a Wilsonian spirit upon integrating out quantum 
fluctuations momentum shell by momentum shell. For this we employ the formalism based 
on the exact renormalization group flow equation for the effective average action 0, 



where denotes the second functional derivative of the effective average action F^ that 
governs the dynamics of the system at a momentum scale k. The logarithmic scale pa- 
rameter t is given hj t = In/c/A, dt = k{d/dk), where A denotes the ultraviolet (UV) 
scale at which we define the bare action Fa. The cutoff function Rk is to some extent 
arbitrary and obeys a few restrictions j3] which ensure that the flow is well defined and 
interpolates between the bare action in the UV and the full quantum effective action F^^q 
in the infrared (IR). 

We solve the flow equation (jS)) by using Eq. (jlj) as a truncation of the space of all 
possible action functionals. As a consequence, we promote all couplings and wave function 
renormalizations occurring in Eq. (0)) to /c-dependent quantities. Although the truncation 
(jD) represents only a small subclass of possible operators generated by quantum fluctua- 
tions, it is able to capture many physical features of QCD-like systems. 

Let us elucidate the single components in detail: for the scalar potential, we use the 
simple truncation 



Already the p-dependent first two terms of the potential are capable of describing spon- 
taneous xSB of the system which we are aiming at. Indeed, the order parameter a de- 
notes the minimum of the scale-dependent effective potential Uk for A; ^ 0. The term 
~ C = + 0* breaks the Ua(1) symmetry of simultaneous axial phase rotations of scalars 
and fermions; it accounts for the effects of the axial anomaly. However, the presence of 
the axial anomaly is not relevant for universality of spontaneous xSB , although it has, 
of course, a strong quantitative impact on resulting low-energy parameters such as con- 
densates and constituent quark masses. Therefore, we postpone the discussion of this 
quantitative influence to Sect. Eland set u = in the following for the sake of clarity. 

In the gauge sector, we do not attempt to calculate the full nonperturbative flow of Zp, 
or alternatively the gauge coupling g, here, but study various possibilities for these flows 



atFfc = ^STr[5ii?fc(Ff 
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(5) 



(6) 
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and take over nonperturbative results from the literature. The most important features of 
the universality classes involve only the perturbative running of g} 

We will define the quantum theories by fixing the initial conditions for the renormal- 
ization flow at the UV scale A. In the gauge and fermion sectors, we choose 

^fL=a = 1' ^■V'L=a = 1' A^L=a = 0. (7) 

The first two conditions normalize the gauge and fermion fields and imply that g denotes 
the bare gauge coupling. The last condition states that four-fermion interactions either 
have been partially bosonized into the scalar sector or are completely absent at the UV 
cutoff scale A. 

The choice of the scalar couplings at the UV cutoff will finally determine whether we 
are in or beyond the QCD domain. In order to describe standard QCD in our picture, a 
natural choice is given by 

m'\,^^ = +OiA'), A^|,^^ = 0, iZ^,h)\,^^^0, (8) 

implying that the scalar fields are nondynamic, noninteracting and heavy at A and decouple 
from the fermion sector. They could be integrated out without any effect on the fermion 
sector and therefore are completely auxiliary. However, we will demonstrate below that 
the infrared physics including xSB is to a large extent independent of the initial values in 
the scalar sector; in other words, the QCD universality class is actually much bigger than 
the restrictive choice of initial conditions of Eq. (jH)).^ 

For a concise presentation of the RG flow equations of the single couplings, it is conve- 
nient to introduce the dimensionless, renormalized and /c-dependent quantities, 

- \ - ^ h - ^ (Q-\ 

in the symmetric regime of the system. In the x^B regime, the mass term becomes negative, 
and we replace this coupling by the minimum of the potential po and its corresponding 
dimensionless variable k defined by 

= |f/.(p = Po), K = ^. (10) 

Similarly, we define as the second p-derivative of the potential at the minimum in 
the xSB regime. The running of the wave function renormalizations is studied using the 
associated anomalous dimensions, 

?70 = -dt In Z^, r]^ = -dt In Z^, r]F = -dt In Zp, (11) 



^Thc running of g is universal in two- loop order. In the framework of the exact renormalization group, 
this has been computed in jHj. As discussed in the introduction, the one-loop running is actually sufRcient 
to generate the main qualitative features needed for our argument. 

^Already at this point, it is clear that A^^a could also be chosen nonzero, which would only result in an 
unimportant change of the normalization of the functional integral. 
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where rjF represents the major piece of information from the gauge sector in our truncation. 
Here, the use of the background-field method for this gauge sector has two advantages: first, 
it represents a book-keeping device to set up consistent gauge-invariant approximations 
within a certain order of truncation. Second, the physical idea of the background field 
is that it accommodates the true ground state of the system around which the quantum 
fluctuations are integrated out. In this spirit, we deduce the running gauge coupling from 
the RG behavior of the background field. Owing to background gauge invariance, the 
product of gauge coupling and background gauge field is renormalization-group invariant 
[Sj, so that the beta function for the renormalized running gauge coupling g is related to 
r]F by 



(3,2 ^dtg'' = VF9\ 9' = ^- (12) 

Actually, the effective action depends on both the background and the fluctuating gauge 
field, and the n-point functions can only be extracted from the functional depending on 
both fields [B]. Nevertheless, once all fluctuations are integrated out, the fluctuating field 
can be set to zero and the resulting effective action is gauge invariant. In general, the 
dependence of the effective action on both fields is needed for the RG flow. With the help 
of background-field identities, the dependence of the effective action on the fluctuating 
gauge field and the background field are related. A detailed record of the flow equations 
and results in the background field formalism, including the role of the gauge symmetry 
and Slavnov- Taylor identities, can be found in [BllZj. 

In the present work, we neglect possible differences between the RG flow for gauge 
couplings defined from the background-field effective action and from vertices of the fluc- 
tuating field jHl E]- This is perfectly justified in the limit of small gauge coupling which 
is of primary importance for this work. Here the lowest-order running is universal. By 
contrast, in the region of large coupling, our truncation of the gauge sector would anyway 
not be reliable if taken at face value, so that we abstain from resolving the gauge field 
running and r]p in this regime. In this region, we simply model the running of the gauge 
coupling in order to obtain a first glance at the xSB regime. Thereby, we assume that the 
influence of higher gluonic operators can be effectively accounted for by the increase of the 
gauge coupling. Although this certainly represents an oversimplification, let us stress that 
the details of the flow in the gauge sector are only of secondary importance for the issue 
addressed in this paper. 

Inserting the truncation into the exact RG flow equation for the effective average 
action, we find the following results. The scalar and fermion anomalous dimensions can be 
written as 

?70 = 4w4KAjm^2(O'2KA^;?70) 
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mf^'^{Kh^]ri^) + nh^ m'2^''^{Kh?']ri^) , (13) 
= 2C2{N^)vi (3 - ^{^^'^(k/i^ 0; r]^, ?7f) - 3(1 - m'ii^'^{Kh^, 0; rj^, 7]^) 
+W4 h"^ [m[™^'^(/t/i^, e + 2KXf, r]^, r]^) + mfp'^^^nh^, e; t]^, , (14) 
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where U4 = l/(327r^) and C2{Nc) = {N^ — l)/(2Nc). This representation is vahd in the 
symmetric as well as in the xSB regime. In the former, k has to be set equal to zero, 
whereas e = has to be chosen in the latter. The various quantities denoted by m 
are threshold functions which control the decoupling of massive modes for decreasing k; 
they also contain all dependencies on the precise choice of the cutoff function i?^. Their 
definitions and explicit representations can be found in App. |X]or in j^. 

Equation (fT^ agrees with |^ and m. We also find agreement for the second line of 
Eq. fll4p . whereas the first line arises from the gauge- field sector (which has not been dealt 
with in |l],jT5). As a further check, we note that in the perturbative small-coupling limit, 
where the threshold functions m occurring above universally reduce to 1, we obtain 

I ^= I ,92m 2^ 1 -2 



which agrees with the literature 

In the symmetric regime, the flow of the purely scalar sector can be summarized by 

dte = -{2-r]^)e-8v4X^lt{e-r]^)+8N,v^h^lf^'\0-r]^), (16) 
dtX^ = 2r]^X^ + 20v^Xll^^{e-r]^)-SN,VihUf'''\0-r]^), (17) 

whereas in the xSB regime, we find 



dtX^ = 2r]^ X^ + 2v4 Xl v<p) + 18^4 l^^nXf, r]^) - SN.v^ l^P'^f^h''; r/^), (19) 

in complete agreement with the results of Again, the quantities denoted by / are 

threshold functions [121 . Now we turn to the flow of the Yukawa coupling, which is 
driven by all sectors of the system: 

dth'^ = (2?7^ + 770) /i^ - 4^4 [/^™^'^(k/i^, e; r/^, rj^) - Z^™^'^(k/i^, e + 2kA<^; r/^, 77^)] 

-8(3 + 0C2(iVc)t'4 g'h^ lfr^\^^h\ 0; r/^, r/p), (20) 

where we have to set k = (e = 0) in the symmetric (xSB) regime. As a check, we take a 
look at the perturbative limit, 

I pert. 87^2 47r2 ^ ' ^ ' 

where we rediscover known results and also observe that the gauge-parameter ,^-dependence 
has dropped out as it should. 

A crucial ingredient is the flow of the fermion self-interaction, which - in dimensionful 
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representation 



can be written as 



dAa 




(22) 



Here we neglected terms ~ k which arise only in the broken regime but are suppressed 
therein owing to simultaneously occurring threshold functions (these terms are similar to 
the last term in square brackets in Eq. (fT^ . which has hardly any effect on the results 
either). In Eq. (j22p as well as in all equations above, we neglected terms of order Ag- 
on the RHS, because = will finally be guaranteed on all scales as discussed below. 
Furthermore, we have chosen the same Fierz transformations in the Dirac algebra as in ^ 
and decomposed the possible color structures of the four-fermion interaction into a color 
singlet (S-P)s and color A''^^ — 1-plets (S-P)Ar2_i, (V)Ar2_i. In the present work, we focus 
on the (S-P)s term; in principle, the (V)7v|-i term could be absorbed into a fc-dependent 
transformation of the nonabelian gauge field in the same way as suggested in j2j for the 
abelian case.^ 

As mentioned above, there is a certain redundancy in the parametrization of the effec- 
tive action owing to possible different choices of partial bosonization of the four-fermion 
interaction. From a different viewpoint, this redundancy corresponds to the possible mixing 
of fields or composite operators with identical quantum numbers. We remove this redun- 
dancy in the present truncation with the aid of the following /c-dependent transformation 
of the scalar field ( "fermion-boson translation"): 



with a priori arbitrary functions <yk{q) and /3fc(g). Upon this transformation, the flow 
equations given above receive additional contributions ~ ak{q),[3k{q) according to 



As described in more detail in j2] , these functions can be uniquely determined by demanding 
for (i) dtXaiq^) to vanish for all k and g^, where the momentum dependence of A^ has been 
studied in the s channel for simphcity, Xaiq"^) = X^{s = g^), (ii) the Yukawa coupling h 

^By neglecting some of the four-fermion interactions, our quantitative result will depend slightly on the 
choice of the Fierz decomposition. Using "fermion-boson translation" to be described in the following, this 
dependence can be removed in a larger truncation, as was recently shown in jlOj . However, we checked 
explicitly that quantitative results in another natural Fierz decomposition involving (S-P)s, (V)s and 
(V)jv2-i differ from the present ones only on the 1% level. 



dt(t)k{.q) 

dMq) 



-(^LV^R)(g) dtak{q) + 0fc(g) dtPkiq), 
(^R^L)(-g) dtakiq) + (pliq) dMq), 



(23) 




(24) 
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to be momentum independent, and (iii) dtZ^i^q"^ = k'^) = —rj^Z^ in order to render the 
approximation of a momentum- independent self-consistent. Condition (i) together with 
the initial condition ((Zj) guarantees that no four-fermion interaction of this type is generated 
under the flow; this interaction is bosonized into the scalar sector at all scales k. Condition 
(ii) guarantees the fermion mass generated by x^B is also momentum independent, so that 
the couplings in the xSB regime have a direct physical interpretation. 

The field transformation (f^Hj) affects also the scalar couplings, and we obtain in the 
symmetric regime: 

d,h' = d,h'\^^ + 2{l + 2e + Q^{l + eY){p(g' + Pllh'), (25) 

where the corresponding first terms on the right-hand sides denote the flow equations for 
fixed fields as given above in Eqs. (fTBj) and ((20)). In the xSB regime, we find similarly 

d,K = + ^ ^^^ ^2^'^ (1 + (1 - ^^^)Q'^) « 9' + h'), 

dth^ = dth\+2{l-2K\^ + Q,{l- KX^f){pfj^ + P^^h"). (26) 

Defining AAo- := A<j(fc^) — Acr(O), the quantity Qa = dt/^\a / dtXaiS^) measures the sup- 
pression of Ao-(s) for large external momenta. Without an explicit computation, we may 
conclude that this suppression implies Q^j < 0, in agreement with unitarity; furthermore, if 
the flow is in the xSB regime, the fermions become massive, and non-pointlike four-fermion 
interactions in the s channel will be suppressed by the inverse fermion mass squared.^ 
Therefore, we model Q„ by the ansatz 

= Ql mS™)'"(ft:/i2, 0, 7]^, rip), Q° = const. < 0, (27) 

where we have introduced a threshold function with the appropriate decoupling properties 
for massive fermions. The qualitative results are independent of the precise choice of Q^, 
and it is reassuring to observe a quantitative independence of the IR observables on the 
precise value for (e.g., ~ —0.1). 

The field transformations (j23p also modify the equation for A^ via the terms ~ dtPk- 
In the pointlike limit (g^ = 0), the modified running is given by 

dtX^ = dA^\^^ + 4^(1 + e)(l + (1 + e)Q^) (/?£ + (3^^ h'). (28) 

It will turn out that the modification of the flow of A,^ is also quantitatively irrelevant, 
whereas the the modifications displayed in Eqs. (j^^ and (j^UI) are of crucial importance. 

"^This can be inferred from the heavy-fermion limit of the two-gluon/scalar-exchange box diagram where 
the internal fermion propagators become pointlike ~ 1 /wf . 
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3 Bound-state fixed point 



The universal features of spontaneous xSB in the QCD domain that will be quantitatively 
analyzed in the next section can be traced back to the occurrence of a fixed point for 
the scalar couplings. This fixed point is infrared attractive as long as the gauge coupling 
remains weak and can be associated with a bound state j2]. 

The fixed-point structure can conveniently be analyzed with the help of the coupling 

^ _ e _ Z|m2 

'kW- ^^^^ 

Since we are interested in the domain of weak gauge coupling, for simplicity we can neglect 
the anomalous dimensions in the following. In this approximation and choosing the gauge 
parameter ^ = (background Landau gauge), the flow of e yields: 

dtl = 8N,v/P^' - SvJtie) ^ - (2 - 2AC,{N,)v/;,^^^' e 

-2{Pllg' + P^lh')i'. (30) 

(Here, all arguments of the threshold functions which are not displayed are assumed to 
be equal to zero; therefore, threshold functions without any argument are simply numbers 
which depend on the details of the regularization). If the scalar field is auxiliary at the 
UV scale as in the QCD context, its wave function renormalization is very small initially, 

-C 1, so that the dimensionless renormalized mass is very large, e ^ 1. In this case, 
scalar fluctuations are suppressed and the threshold functions depending on e vanish; the 
right-hand side of Eq. (j3(jp describes a parabola in the variable e, and we find two positive 
fixed points, < < ^2, where is UV attractive but IR unstable, and is an IR stable 
fixed point (see Fig. |21 solid line). It can be shown that el corresponds to the inverse of the 
critical coupling of the NJL model, so that our flow describes a model with strong four- 
fermion interaction if we choose UV initial conditions with eA < el to the left of the first 
fixed point (see, e.g, for a detailed analysis of the phase structure in the abelian case). 
For this choice, the system is not in the QCD domain but approaches chiral symmetry 
breaking (e < 0) in a perturbatively accessible way (P^SB). 

In this section, we concentrate on those initial values which release the system to the 
right of the first fixed point, eA > Ci, i.e., which are weakly coupled in the NJL language. 
This will be the range of the QCD universality class. As the system evolves, it flows 
towards the second fixed point which then governs the evolution over many scales. 
Here, the system "loses its memory" of the initial conditions; in particular, it is of no 
relevance whether we start with < eA < e2 or ca > 62- The evolution towards and in 
the IR is universally governed by this fixed point eg, which can be shown to be associated 
with a fermion-antifermion bound state; e.g., in QED, the properties of the scalar field at 
this fixed point correspond to those of positronium [2]. 

Before we elucidate the fixed-point properties further, let us briefly mention that its 
existence can be generalized to the case of a scalar field describing a fundamental particle 
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Figure 2: Flow of e according to Eq. ()30|) (schematic plot): the solid line corresponds to a 
QCD scenario at weak gauge coupling; the arrows indicate the direction of the flow towards 
the infrared. The dashed line corresponds to a system with fundamental scalar, Z^\k=/^ = 1, 
e < 1, and strong scalar self-interaction. The dotted lines exhibit the destabilization of the 
bound-state fixed point by the increasing gauge coupling. 

in the UV (a Yukawa model with gauged fermions rather than QCD). In this case, we 
have = 1 and e ~ 0{1) at the UV scale. Now the second term in Eq. (jHUj) can become 
important, in particular for a large 0^ coupling and/or small h"^. When discussing the 
RHS of Eq. (jHUj) for fixed g, h, one should keep in mind that these couplings may 
change with k. For large X^/h"^, the e parabola is lowered and the first fixed point can 
move to negative values, < (see Fig. |2l dashed line). In this case, we can release the 
system even in the broken regime at the UV scale, e, e < 0, but it still evolves towards the 
bound-state fixed point eg. In comparison with Fig. Q this corresponds to initial values 
< fh\ < 0. Physically, such a scenario describes a system involving fundamental 
scalars, fermions and gauge fields, where the scalar sector is initially weakly coupled to 
the fermions. If we start in the broken regime, scalar fluctuations will drive the system 
towards the symmetric regime before the fermion-gauge-field interactions induce sizable 
bound-state effects which can exert an influence on the scalar sector. In this scenario, the 
first fixed point < is a measure of the strength of the initial effective coupling between 
scalars and fermions. For strong effective coupling, < el, an initial negative scalar mass 
of the order of the cutoff, ~ — O(A^) will induce a vacuum expectation value and 

a fermion mass of the same order, in agreement with naive expectations. But at weak 
effective coupling, e.g., h'^ ~ C'(l), A,^ ~ 100 and < eA < 0, the system can still start 
with an initial negative scalar mass rh'^\j^_^ — — (9(A^), but finally run into the bound-state 
fixed point. As an important result, the vacuum expectation value and the fermion mass 
after symmetry breaking can easily be orders of magnitude smaller than the UV scale, as 
exhibited in Fig. ([Q) in the Introduction. We conclude that all systems with eA > belong 
to the QCD universality class. 

Let us now turn to the properties of the system at the bound-state fixed point. The 
crucial observation is that not only e but also all dimensionless scalar couplings approach 
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fixed points. In the general case, the fixed-point values depend in a complicated form on 
all parameters of the system. However, in the limit e ^ 1 (QCD-like), we can find analytic 
expressions that satisfy the fixed-point conditions dt{e, /i^, A^) = to leading order: 



2 

5 



~ ^'^Cl^' _ 12 C2(iVc)(iVc + 2)(iVe-l) ,(FB),4 4 



(31) 



A* 



From the first equation, we read off that the approximation e ^ 1 is equivalent to assuming 
\Qa\ ^ 1, which is roughly fulfilled in our numerical study with our choice of Q° = —0.1. 

The remarkable properties of the IR fixed point become apparent when considering the 
renormalized scalar mass, = e/c^. Since e — *• e*, the scalar mass simply decreases with 
the scale k, so that it is only natural to obtain small masses <C rh\ for small scale 
ratios A; <^ A. In other words, even if we start with a scalar mass of the order of the cutoff, 
'^^l/j-A ~ iio fine-tuning will be necessary to obtain small mass values at low-energy 
scales, as long as the running is controlled by the bound-state fixed point. 

In order to approach the xSB regime, the bound-state fixed point has to be destabilized; 
otherwise, the system will remain in the symmetric regime as is the case in QED. In QCD, 
this destabilization arises from the increase of the gauge coupling towards the infrared J3] . 
From the third and last term of Eq. (|3(jp. it is obvious that an increasing gauge coupling 
lifts the dfi parabola (see Fig. |21 dotted lines). For some value of the gauge coupling, 
the two fixed-points in e will be degenerate, so that there is no fixed point at all for all 

> g^. The beta function dte is then strictly positive, which drives the system towards 
the xSB regime. 

In the limit e ^ 1, the critical gauge coupling of fixed-point degeneracy g^ can be 
computed analytically, and we find: 



P 16 , iVc / / 1 \ 4 . 1 

where we have used linear cutoff functions ^Hj for which = 3/2. For instance, for 

SU(3) we get q;d = 4^ — f, which is in the nonperturbative domain, as expected. As soon 
as g^ exceeds g^, the running of the scalar couplings is no longer protected by the bound- 
state fixed point. Here all couplings are expected to run fast, being strongly infiuenced 
by the details of the increase of the gauge coupling. Of course, owing to strong coupling, 
many higher-order operators can acquire large anomalous dimensions and contribute to 



■^Strictly speaking, the value of is not a physical quantity and depends on the choice of the cutoff 
function, i.e., regularization scheme. This is only natural, since the running of the coupling itself also 
depends on the regularization. The scheme dependence, however, cancels out in physical quantities. 
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the dynamics of the symmetry-breaking transition. Our truncation should be understood 
as the minimal lowest-order approximation in this regime, but gives already a remarkably 
consistent (but not necessarily complete) picture. Once chiral symmetry is broken, the 
fermions decouple and the fermionic and (most of the) scalar flow essentially stops. 

The scenario discussed here finally explains why the IR values of the scalar and fermionic 
couplings inherit their order of magnitude from the QCD scale Aqcd as they should, 
whereas particularly the details of the scalar sector at the UV scale are of no relevance, 
owing to the fixed-point structure inducing QCD universality. 



4 Numerical results 

In the following, we concentrate on the set of theories that belong to the QCD universality 
class. In order to illustrate how universality arises from the presence of the bound-state 
fixed point, we initiate our flows at a GUT-like scale of A = lO^^GeV, where the gauge 
coupling is weak and increases only logarithmically towards the infrared. Therefore, the 
bound-state fixed point exists over a wide range of scales. As discussed before, hardly any 
dependence on the specific initial values for the scalar potential and the Yukawa coupling 
remains because of the fixed point, as we will demonstrate quantitatively in the following. 

For illustrative purposes, we concentrate here on QCD-like scenarios where the scalar is 
auxiliary at the UV scale, and explore this parameter space using the natural choice given 
by Eq. (jHl) as a reference; to be precise, we use the reference set, 

,21 _ a2 \ I _ n V I _ in-8 121 _ i n-12 



Ia 



^ el =10^ A^|a = 0, /^|A = 10-^ (33) 



in our numerical studies. In all computations, we use linear cutoff functions proposed in 
[TB] for which the threshold functions can be determined analytically (see App. We 
plot the flows of the renormalized dimensionless couplings e, h and A,^ in Fig. |21 for the 
symmetric regime. The reference set ()33|) is depicted as solid lines, whereas the dashed 
and dotted lines correspond to initial values which deviate from the reference set fj33|) by 
many orders of magnitude for the corresponding couplings. 

As long as we start in the range of attraction of the bound-state fixed point, we can 
obviously vary the initial values for the scalar couplings over many orders of magnitude 
without any appreciable effect. The system quickly approaches the bound-state fixed 
point, where the initial values of the couplings become unimportant. In particular, the 
scalar mass, which is allowed to be of the order of the cutoff or even much larger at /c = A, 
runs to small values ~ k while the system is governed by the bound-state fixed point. No 
fine-tuning is necessary for this.^ Let us stress once more that these features of universality 
are not restricted to the reference set ()33p and the variations thereof. They can also be 

^As a fairly weak condition, we only have to ensure that the initial scalar couplings are such that no 
strong four-fcrmion interaction is implicitly induced by the initial values; in this case, the system starts to 
the left of the first fixed point and rather resembles a nonabelian gauged NJL model with PxSB . 
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k = 10*10 GeV 

Figure 3: Flow of e, h and A,^ in the symmetric regime according to Eqs. (|16p. (|17|). 
(l2Up . and (j25|) . The sohd hnes correspond to the reference set (j33p . whereas the dotted 
and dashed hnes represent the flows for strongly differing initial values as indicated. The 
insensitivity with respect to the choice of initial conditions is clearly visible. On the 
horizontal axis, the exponent tio is used for the scale k = 10*"^GeV. 

found in Yukawa models with a fundamental scalar {Z^\k=\ = 1) and even if we start in 
the broken regime at the U V scale (see Fig. . 

At the bound-state fixed point, the couplings are modulated only by the logarithmically 
slow increase of the gauge coupling. Incidentally, the modulation of e = e/h'^ is completely 
carried by h, whereas e stays fixed. This agrees with our analytical fixed-point values found 
in Eq. (j31|) . A rapid change for the couplings in Fig. |21is visible after g"^ exceeds and 
the bound-state fixed point has disappeared {tio < 1). 

The behavior of the system changes rapidly after the gauge coupling has grown large. 
For g"^ > g'^, the bound-state fixed point vanishes and all couplings start to run fast. 
The system necessarily runs into the xSB regime where the scalars develop a vacuum 
expectation value and the fermions acquire a mass 

= lim P kK^ = {ha^f, (34) 

where ctr = limfc^o \/Z<i> Po denotes the renormalized expectation value of the scalar field. 

This leads to a decoupling of the fermions, and, consequently, fermion-boson translation 
is "switched off". Also the flow of the Yukawa coupling stops, the scalar and fermion 
anomalous dimensions approach zero, and k runs according to its trivial mass scaling, 
K ~ 1/fc^, so that mf approaches a constant value. 

Whereas the qualitative picture is rather independent of the details of the running gauge 
coupling, quantitative results are highly sensitive to the flow of the gauge sector. This is 
because a finite amount of "RG time" passes from the disappearance of the bound-state 
fixed point to the transition into the x^B regime. In between, the running of the gauge 
coupling exerts a strong influence on all other couplings which are no longer protected 
by any fixed point. A purely perturbative running of the gauge coupling turns out to be 
insufficient for the present purpose, since the (unphysical) Landau pole destabilizes the 
system in the infrared. 
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For definiteness, let us consider a running coupling governed by the beta function 



for our numerical studies. In the UV, this beta function exhibits an accurate two-loop 
perturbative behavior, whereas the coupling runs to a fixed point as = 5'^/(47r) — a^, in 
the IR for k —>■ 0. In the first place, the infrared fixed point is convenient for numerical 
purposes, since it does not lead to artificial IR instabilities. Moreover, an infrared fixed 
point for a mass-scale-dependent running coupling is compatible with the expected mass 
gap in Yang-Mills theory. Below this mass gap, all gauge field fluctuations decouple from 
the flow and can no longer drive the flow of the coupling. Different beta functions with 
and without infrared fixed points are studied in Appendix IHl It turns out that, though 
the infrared properties such as the constituent quark mass depend quantitatively on the 
choice of the beta function as expected, the universal features discussed in the following 
remain untouched. This underlines our observation that the detailed understanding of the 
flow for the region of strong gauge coupling is not essential for the overall picture. 

In combination with Eq. (|35|) . the system of flow equations is now closed and provides us 
with an answer for the (truncated) quantum effective action, once we specify all parameters 
and initial values. We have investigated SU(A'c = 3) gauge theory with initial value g{A) 
chosen such that acquires its physical value at the Z-boson mass, as{Mz) — 0.117. We 
work in the background Landau gauge, ^ = 0, which is known to be a fixed point of the 
renormalization flow in the gauge sector P7j,^H]- If we had an exact flow equation at our 
disposal this choice would fix the system completely. 

In our truncation, however, we have the parameter in addition to the Yang-Mills 
beta function, which characterizes our ignorance of the exact flow. The quantity mea- 
suring the momentum suppression of the four-fermion interaction will be set to = —0.1, 
in agreement with our considerations given above. It turns out that the infrared properties 
of the system are only weakly dependent on this parameter and on (see below), which 
substantiates our truncation. Furthermore, we choose a* to be of order 1, but not too close 
to g^/{4:iT) in order to avoid pathologies: a* = 2.5. 

For this concrete scenario, the transition to the xSB regime occurs at k^sB — 423MeV. 
The renormalized scalar mass slightly above k^sB and the VEV of the scalar field below 
k^SB are depicted in Fig. El (left panel). According to Eq. we find a constituent quark 
mass of TTif ~ 371MeV as shown in Fig. |3] (right panel). Of course, these numbers depend 
strongly on the details of the Yang-Mills beta function for strong coupling ~ 1; various 
other examples are discussed in Appendix El Finally, the running of A,^, h"^ and the scalar 
and fermionic wave function renormalizations is collected in Fig. El 

Focusing on low-energy QCD-like aspects of our truncated system, it is also remarkable 
that (apart form the scalar couplings) the choice of has little effect on infrared properties 
of the system: varying between —0.5 . . . 0.001 changes k^sB or rrif only at the level of 
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Figure 4: Flow of the scalar mass m, the scalar VEV ctr, and the constituent quark mass 
rrif close to and in the xSB regime, using the reference set flHHj) . For the particular choice 
for the running of the gauge coupling according to Eq. with a^, = 2.5, the transition 
occurs at fc^^sB — 423MeV. 



less than 10%. This is reassuring and in contrast to the strong Q^-dependence of the 
bound-state fixed-point values of and The variations of the infrared properties are 
similarly small for changes in the gauge parameter in the interval ^ = ... 2. 

To summarize, a large class of Q CD-like theories including a scalar degree of freedom 
belong to the QCD universality class owing to an attractive infrared fixed point present 
for weak gauge coupling. Even before the gauge coupling becomes strong, all theories 
in this universality class are indistinguishable at low energies. They exhibit an identical 
approach to x^B which is triggered and quantitatively determined by the increase of the 
gauge coupling. 

5 Instanton-mediated interactions, axial anomaly 
and the fate of the eta boson 

Up to now, we have considered only that part of the model which has a global Ua(1) 
symmetry corresponding to simultaneous axial phase rotations of the scalars and fermions. 
In QCD, this symmetry is anomalously broken by the presence of gauge-field configurations 
of nontrivial topology. For instance, instantons induce fermion interactions which break 
this symmetry. In an instanton-anti-instanton background, the A^f = 1 interaction is mass- 
like and can be expressed as [12] 

A+A = / 4 do'iQ) CE(iVc) (27r2^,3) JJJ^ ^^^^^ _ ^^^^)^ (36) 



where Ce{Nc) is a color factor that arises from averaging over all possible embeddings of 
SU{2) into SU{Nc), e.g., Ce(2) = 1, Ce(3) = 2/3, and /i=lGeV is the renormalization 
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Figure 5: Flow of A^, /i^, and the wave function renormalizations and over the 
complete range of scales for the reference set ()33|). The rapid change of all couplings near 
tiQ = \ogiQ k^sB /A ~ —0.5 is visible. Whereas h"^, and Z^ approach fixed points in the 
deep infrared owing to decoupling, decreases logarithmically owing to a massless "eta" 
in absence of the axial anomaly. 



scale for the fermion fields. Note that we introduced an IR cutoff function fc{k,m{) in 
the upper bound of the instanton radius g integration. This function should cut off the 
contribution from all modes with momenta either below k or the generated fermion mass 
mf, and thereby implements the renormalization group formulation of this interaction in a 
simple manner. The g integration is UV finite for g ^ owing to asymptotic freedom, and 
the infrared {g oo) is controlled by the cutoff fc and by the increase of the coupling. 

In the following, we intend to include this interaction as it is, being an example for 
a Ua(1) violating term. Contrary to standard instanton based models [20]; we do not 
employ further information about, e.g., average instanton sizes and separations or other 
assumptions about the vacuum state of the gauge field. For this, we note that Eq. 
already corresponds to an integrated flow, £i+a = ?7ii+A (^r^l — ^/'l^/'r), where the flow of 
the induced mass mj+A is given by'' 



dtfcik^rrii), (37) 

Q=Ic{k,mi) 

more rigorous treatment of anomalous Ua(1) breaking within the flow equation formalism has been 
suggested in |21| . 
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with the initial condition rni_|_A(fc = A ^ oo) 0. For consistency, we also included 
here the fermion wave function renormalization Z^, which was not taken into account in 
Eq. p6|) as derived in ^H]- Since fc has mass dimension -1, an appropriate choice is given 
by 

Uk,m,) = lfciKh'), with /,(0) = 1, (38) 

such that /c(0,mf) = l/rrif. For our numerical solutions, we will use fc{x) = (1 + x)^^^'^ 
for simplicity. With these definitions, we can rewrite Eq. (jHTj) as 

dmi+A = -2n'Z^ A d^^ifjk) CE(iVc) ( ^it/^) i 1 + d,{^h^)\ , (39) 



/c ^ V— ---- ^ «(^) J y 

where fc = fdufi^), and the prime denotes a derivative. 

Now we could repeat the calculation of the flow equations of SectO including this 
fermion mass term in the propagator. In this way, however, we would induce a number of 
Ua(1) noninvariant fermion- fermion and fermion-scalar couplings which complicate the cal- 
culation unnecessarily. Instead, we propose a generalization of the field transformation 
which serves to translate the instanton-induced interaction into the scalar sector: 

dt4>k{q) = -{i^Li^R){q) dtak{q) + <pk{q) dt(3k{q) + d^k + {(p*k4>k)(f>kdtSk, 

dtKiq) = (V^RV^L)(-g)5iafc(g)+0K9)5t/5fc(g) + %;^ + (</':</'fc)0:5A, (40) 

with additional a priori arbitrary functions 7^ and 5ki whereas Uk and 13k are those of 
Sect. 121 The term ~ dt'^k corresponds to a Ua(1) violating shift of the scalar field which 
can compensate for the instanton-induced fermion mass. The flow of mi+A is now given 
by ^ 

dtm+A = dtThi+A\^^ + h dak - -i" dtak, (41) 

where the second and third terms arise from the transformation of the Yukawa interaction 
and the last term in Eq. (jU)), respectively. Now we can determine 7^ such that dtfhi^x = 
holds on all scales. In this way, the instanton interaction does not affect mj+A (which 
vanishes on all scales), but is translated into the scalar sector and contributes to the 
running of v. In the point-like limit (g^ = 0), we find: 

dtv = -2m^dak + ydtf3k. (42) 

Introducing the dimensionless renormalized quantity 

.= ^, ^ u^ = k^.^ (43) 
where i^r denotes the renormalized (dimensionful) value, we finally arrive at 

+^(l + (l + 6)U)(/5f>^ + /3f>^), (44) 
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which describes the running of the axial anomaly in the instanton approximation. 

The shift ~ dt'^k induces another Ua(1) violating term {(jf (j)){(f)* + 0) via the transfor- 
mation of the X^{(j)*(j))'^ term. This can be cancelled by an appropriate choice of the last 
transformation function 6k in Eq. (jlOI), which has to satisfy 

A</,<9f7fc - = 0. (45) 

Finally, the terms ~ Sk in Eq. (jlHI) influence the running of A<^ via the transformation of 
the scalar mass term. The modified flow equation for A<^ reads: 

dtX^ = dtX4^^+A^{l + 2e+il + efQ,){plg' + f3^lh') 

+167r2^rf^^(/,A)CE(iVe) ( ^^y'^"" ^(l+^-^d,{Kh')] . (46) 



These equations are valid in the symmetric regime with similar equations for the xSB 
regime displayed in appendix O 

Strictly speaking, the system is never in the symmetric regime, since chiral symmetry 
is always broken implicitly by a nonzero u term which induces a nonzero VEV cxo for 
the scalar field. For instance, rotating the VEV into the real component, (p = ao = (p* , 
c"o = a/PO' location of the minimum obeys 

= U'ipo) =m^ + Vo - ^ = e + kA^ - (47) 

^VPo ^ V 



Obviously, k = is not allowed if z/ 7^ 0, owing to the linear term in (j) in Eq. The 
running of the minimum can be inferred from 

= dtU'ipo)\^ = U"ipo)dtpo + dtU'ipo)\^^ dtpo = -^ 7-371 ^*^'(/^o)|,„. (48) 

Since the instanton-induced terms are exponentially small for the major part of the flow, the 
minimum of the potential is actually very close to zero, and the equations for the symmetric 
regime of Sect. |21can be used up to tiny corrections. The solution of the flow equations 
is numerically difficult with an exponentially small k, in the broken regime. Therefore, we 
decide to solve the flow equations for large enough k in the symmetric-regime formulation. 
In this regime, u evolves according to Eq. (jl^ with only a subdominant coupling to the 
other flow equations via Eq. pUj) . Then we switch to the broken- regime description at 
that scale where the instanton-induced fermion mass mf is of the order of a few MeV; this 
procedure induces an error only at the per-mille level and turns out to be insensitive to 
the details of the switching scale. 

We have analyzed the flow equations including the instanton-mediated interaction nu- 
merically and used the reference set of initial conditions as defined in Sect. |3] (see Eq. (I33|) ) 
for a direct comparison. As expected, most properties of the system are unaffected by 
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Figure 6: Left panel: axial anomaly uji in the vicinity of the scale of fermion decoupling. 
Right panel: instanton-induced fermion mass. Both plots refer to the reference set flHH|l 
and the particular choice for the running of the gauge coupling according to Eq. ()35|) with 
a^: = 2.5. Comparison with Fig. |3] shows that the fermion mass is dominated by instanton 
effects. 



the instantons, while the system is governed by the bound-state fixed point. Here the 
instanton-induced effects are exponentially suppressed, since the coupling is small. In par- 
ticular, the running of the scalar mass e and the Yukawa coupling are identical to the ones 
displayed in Fig. |21 and the universality properties discussed in Sect. El remain unaffected. 

The renormalized axial anomaly z/r is plotted in Fig. |S1 (left panel). It remains exponen- 
tially small for a large part of the flow and becomes of order (GeV)^ and larger only in the 
strong-gauge-coupling regime. Here, however, it contributes strongly to the VEV of the 
scalar field and consequently to the constituent quark mass which leads to the decoupling 
of the fermions. 

We observe a rather smooth onset of fermion-mass generation. Furthermore, the con- 
stituent quark mass is strongly enhanced by the instanton interactions. For the reference 
set, we find mf = 1765MeV in the infrared limit k ^ 0. Again, this number depends 
strongly on the precise choice of the running gauge coupling in the infrared, and a number 
of other possibilities including instanton effects is listed in Appendix ^ 

Let us finally discuss the fate of the "would-be" Goldstone boson, which we may call 
the eta boson in the style of real QCD. Neglecting the axial anomaly, this boson arises from 
spontaneous breakdown of the global Ua(1) as a true massless Goldstone boson; its effects 
on the scalar sector even after xSB are visible in the logarithmic running of the scalar 0^ 
coupling as can be seen in Fig. El The Ua(1) anomaly, however, generates a mass of 
the eta boson. In the present formulation, the Ua(1) anomaly occurs as the u term in the 
scalar potential ©• Its contribution to the renormalized eta mass can be computed as 

Within the above-given framework of instanton-mediated interaction, we find for the eta 
boson mass in the QCD universality class a value of ^ 4440MeV. Of course, this value 
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also strongly depends on the choice of the running of the gauge coupling and should be used 
only for comparison with other masses computed for the same running gauge coupling. In 
particular, we find roughly the ratio m^/mf ~ 3. This scenario giving rise to a heavy mass 
of a would-be Goldstone boson is familiar from three-flavor QCD. 

By contrast, the fate of the eta boson is more spectacular if we go beyond the border 
of the QCD domain to that of P^SB, corresponding to a choice of rh\ < fhl in Fig. ^ 
or eA < €\ in Fig. |21 Here, the VEV of the scalar field is generically of the order of the 
cutoff A = lO^^GeV. At the same time, the fermions rapidly become massive and decouple 
from the flow only a little below A. As a consequence, instanton contributions or other 
long-distance topological properties have little effect on the fermion sector and thus the 
axial anomaly exerts hardly any influence on the scalars. As a result, the contributions to 
the eta mass are strongly suppressed - powerlike in the denominator and exponentially in 
the numerator. For instance, for the set of initial parameters corresponding to Fig. ^ (right 
panel) with ffi\ slightly below m^, we find an extremely small eta mass, ~ 2 ■ 10~^°eV. 
For smaller fh\, the eta mass decreases even further, and larger eta masses require a 
tremendous fine-tuning of ffi\ close to m^. 

In this scenario beyond the QCD universality class, we have thus found a mechanism to 
generate extremely small masses without any fine-tuning. From another perspective, this 
mechanism exploits the fundamentally different RG properties of scalars and chiral gauge 
theories. For systems in the universality class of P^SB , the x^B scale of the scalar sector 
is generically of the order of the UV scale, whereas the nonperturbative scale of the gauge 
sector can be much smaller. Now the mass of the would-be Goldstone boson is generated 
by the nonperturbative sector of the gauge theory which is exponentially suppressed at the 
UV scale. This interplay finally leads to the generation of the extremely small mass. 

6 Conclusions 

In this work, we studied a class of theories involving one-flavor massless QCD and a chiral 
color-singlet scalar field. Our model is parametrized by the gauge coupling and a number 
of scalar couplings. In this framework, we identified the QCD universality class of theories 
which share the same physics at low energies, namely spontaneous breaking of chiral sym- 
metry triggered by the strongly interacting gauge sector at the QCD scale. As a remarkable 
result, the QCD universality class contains theories with fundamental scalars where the 
microscopic scalar potential has its minimum at nonzero field {fh\ > ffil ~ —0{h?)). For 
these the theories, the scalar fluctuations drive the system first into the symmetric regime 
with a large positive scalar mass, and the remaining flow is governed by the QCD sec- 
tor. We checked explicitly that this is in accord with perturbative expectations for weak 
couplings (cf. Fig. m right panel). 

The mechanism that establishes QCD universality is the occurrence of an infrared at- 
tractive bound-state fixed point in the scalar couplings which persists over a wide range of 
scales as long as the gauge coupling is weak. At this fixed point, the scalar field exhibits 
quark- ant iquark bound-state behavior and the RG running of the scalar couplings is gov- 
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erned by the RG behavior of QCD. All memory of the scalar initial conditions is lost by 
the system. As a remarkable consequence, the scalar mass is not a relevant operator at this 
fixed point. For increasing gauge coupling, the bound-state fixed point is destabilized and 
the system runs towards the xSB regime. Here the role of the scalar field changes and it 
can characterize (quark) condensates and (mesonic) excitations on top of the condensate. 
At strong coupling, the simple overall picture of xSB arising from our truncation can, of 
course, be modified quantitatively as well as qualitatively by the influence of higher-order 
operators. In particular, mixed non-minimal fermion-gluon and scalar-gluon operators 
might add new features to xSB by providing a coupling to the nontrivial gluonic vacuum 
structure. 

Beyond the QCD universality class, we find the class of theories exhibiting perturbative 
spontaneous chiral symmetry breaking (P^SB). In this class, the system is mainly driven 
by the scalar sector, and IR properties such as condensates and generated fermion masses 
depend strongly on the initial scalar parameters. The gauge sector exerts hardly any 
influence on the fermions in this class unless the scalar parameters are fine-tuned to a high 
precision. In the deep IR, pure gluodynamics without dynamical quarks remains. The flow 
of the scalar couplings is never in the attractive domain of the bound-state fixed point, but 
is governed by a fundamental-particle fixed point. Small deviations from this fixed point 
have an infrared unstable component which corresponds to the RG relevant scalar-mass 
operator. 

In both universality classes, we found interesting implications. Our setup of the QCD 
universality class admits a resolution of an old puzzle: whereas QCD has no fine-tuning 
problem and is completely determined by fixing the coupling at a certain scale, low-energy 
QCD models based on NJL-type fermion self-interactions depend strongly on additional pa- 
rameters such as an intrinsic UV cutoff. In the context of partial bosonization, this cutoff- 
dependence corresponds to a strong dependence of IR observables on the bosonization scale 
(or the value of the scalar mass at this scale). In our approach with scale-dependent field 
transformations, partial bosonization occurs at all scales, and no artificial dependence on 
unphysical scales is introduced. In our truncation, QCD flows continuously from a high 
scale with quarks and gluons as the relevant degrees of freedom to intermediate scales with 
quarks, gluons and quark bound states and further to low scales with constituent quarks, 
condensates and mesons. 

In the PxSB universality class with one fermion flavor, we identified a natural mecha- 
nism for the generation of extremely small scalar masses without fine-tuning. The mech- 
anism exploits the fact that the spontaneous breaking of the Ua(1) symmetry would lead 
to an exactly massless Goldstone boson in absence of the gauge interactions. The axial 
anomaly in the gauge sector then endows a small mass to this boson. Owing to the highly 
different RG behavior of the scalar and the gauge sector, the scale of P^SB differs gener- 
ically from the scale of nonperturbative gauge effects by many orders of magnitude. This 
leads to an exponential suppression of the influence of the axial anomaly and thus to an 
exponentially small but nonzero scalar mass. 

For theories with a fundamental scalar, the question arises as to whether our tech- 
nique of fermion-boson translation is capable of describing all possible mesonic degrees of 
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freedom. Let us first look at two extreme situations. For a large negative renormalized 
scalar mass term, perturbation theory applies: there is a fundamental scalar, and sepa- 
rately propagating meson states may not exist - similar to a very heavy top quark. For 
a positive renormalized mass term, the fundamental scalar decouples from the low-energy 
sector in perturbation theory. The low-energy sector then is QCD without scalars, as in 
our picture. The transition is less obvious: in the region where the fundamental scalar 
mass would perturbatively be of the order of the strong interaction scale, there is a strong 
mixing between operators corresponding to fundamental and composite scalars. 

In principle, in a situation with mixing the propagator in the scalar sector may have 
one or several pole-like structures that can be associated with particle excitations. Our 
truncation cannot fully resolve this issue, since, by construction, it follows the flow of only 
one pole in the propagator. Our investigation shows the consistency of a picture with 
only one pole. If the true physical situation had two poles our truncation would follow 
the flow of the lowest mass. We see, however, no indication that a second pole actually 
exists. Nevertheless, it seems worthwhile to discuss the possible implications of a second 
scalar "pole" for the issue of universality classes. First of all, for m\ larger than but not 
in the immediate vicinity of m^, a "second pole" could correspond only to an additional 
heavy scalar particle. This would decay at a high rate into the QCD mesons, since no 
quantum numbers forbid such a decay. (One expects at most a resonance rather than a 
true pole.) Furthermore, effects from the exchange of such a heavy scalar resonance would 
be suppressed by inverse powers of the mass and therefore play no role for the low-energy 
theory. This is what one usually understands by "QCD universality class" , (a notion that 
is not thought to resolve the detailed short-distance physics). 

This issue becomes more interesting when m\ is fine-tuned to the immediated vicinity 
of m^. In this case, we approach the boundary oi the QCD universality class. We emphasize 
that this boundary is not uniquely defined in terms of the symmetries and particles char- 
acteristic for the QCD universality class. Considering the QCD universality class from the 
viewpoint of a larger space of models or parameters, the spectrum of excitations that are 
relevant at the boundary can depend on the direction in parameter space from which the 
QCD universality class is approached. Different directions may yields a differend "number 
of poles" in the boundary region. For this reason, a future more detailed investigation of 
this issue would be quite interesting. 

We stress that all of our main conclusions can be drawn from a mere perturbative knowl- 
edge of the gauge sector which is well under control. In a broader sense, the pure QCD 
sector in our work can be regarded as a particular example for possible other (nonpertur- 
batively) renormalizable theories leading to fermionic self-interactions in scalar channels. 

Let us finally discuss our findings from a different perspective, concentrating on the 
scalar sector. Scalar fields are known to lead to profound problems in quantum field theory 
for two reasons: triviality and (un-)naturalness. Triviality tells us that an interacting scalar 
theory requires a UV cutoff which cannot be removed without switching off the interaction. 
Therefore, whenever we see a scalar quantum field at some low scale, we know that there 
must be new physics at a higher scale. The problem of naturalness tells us that it is 
difficult to achieve a large separation of scales for models with interacting scalar fields 
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without fine-tuning. 

Our formulation lias the potential to solve both problems. A first example can be 
given within the QCD universality class. Although from a QCD perspective, the scalar 
field could be regarded as purely auxiliary, nothing prevents us from considering it as 
fundamental, since the concepts of compositeness and fundament ality are interchangeable 
from the viewpoint of our flow equation with field transformations. We showed in detail 
that "standard" QCD at low energies is indistinguishable from QCD with a fundamental 
scalar, as long as the latter system is in the QCD universality class. In this way, we can 
circumvent triviality by starting in the UV from a scalar field theory without self-interaction 
and Yukawa coupling for which the continuum limit can be taken trivially. The scalar 
interactions are induced by quantum fluctuations. In this construction, the system is always 
in the QCD universality class, and therefore inherits the number of relevant and marginal 
operators form QCD. In particular, the scalar mass term is not a relevant operator, so that 
no naturalness or fine-tuning problems arise in and from the scalar sector. Alternatively, 
we could also follow the bound-state fixed point to /c — oo, where it presumably becomes 
an exact fixed point even beyond our truncation. (The /5 function for the running gauge 
coupling vanishes, owing to asymptotic freedom in this limit.) 

Perhaps more interesting is a second possibility in the PxSB class. Let us consider for 
/c — ^ oo a scalar model with ^ 0, ^ and m^, h chosen such that e corresponds 
to the fixed point e*. This model has an alternative interpretation as a model with four- 
fermion interactions (and no scalar field). Both the gauge coupling and the critical four 
fermion coupling 

vanish for k oo. If this fixed point persists beyond our truncation, it defines a nonper- 
turbatively renormalizable theory ]22^- For lower k a nonzero is generated by the flow 
and we end up with a theory that effectively looks like a model with an interacting funda- 
mental scalar field. This scalar field can give mass to the quarks by P^SB independently 
of the strong interactions, in analogy to the Higgs scalar. The triviality problem could be 
solved in this case - but not the naturalness problem, since we expect a relevant parameter 
corresponding to the scalar mass term. 

This discussion sheds new light on the continuous transition between the P^SB and 
QCD universality classes. In the language of statistical physics, it can be considered as 
a type of crossover between the "fundamental fixed point" and the "bound-state fixed 
point" €2- As a particularity, the gauge coupling is a marginal parameter for both fixed 
points. The scale where it becomes strong sets the lowest possible scale for the effective 
fermion masses. 

Quite generally, the existence of a bound-state-like fixed point leads to a mechanism 
with a naturally small scalar mass. In a sense, this is a realization of earlier ideas of a 
large anomalous mass dimension for the scalar field or "self-organized criticality" [22] • It 
would be interesting to know if a similar mechanism could contribute to an understanding 
of electroweak symmetry breaking which occurs at a characteristic scale hundreds of times 
bigger as compared to QCD. 
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Appendix 



A Threshold functions 

The regularization scheme dependence induced by the cutoff function i?^ is carried by the 
threshold functions / and m. Let us represent the cutoff functions in the scalar, fermion 
and gauge sector by 



q- 

(A.l) 



where y = q^/k"^, and r and rp denote dimensionless cutoff shape functions. Furthermore, 
it is useful to introduce the inverse average propagators P{x) = x{l+r{x/k'^)) and Pf{x) = 
x{l +rp{x/k^))'^, where x = q^. 

Most of the threshold functions given above are defined in Appendix A of |3]. The ones 
which cannot be found therein are marked with a tilde. These can be defined as follows: 



~ (FB),a/ X 



--k'-'' r dxx'/^-^~dt 
2 Jo 



1 + rp{x/k'^ 



2 



-'^ ' dxx'^/^-^dt 

'0 



r{FBB),a/ X 

'1,1,1 («^F, wbi, wb2] rj^, Vb) 



2 



dxx'^/^-^d, 



Pf{x) + k'^wp P{x) + k^Wi 

Pf{x) 1 

(Pf(x) + k^wp)^ (P(x) + k^WB^ 

Pf{x) 1 1 

{Pf{x) + k'^WFy P{x) + Fwbi P{x) + k'^WB2 



(A.2) 



(A.3) 



, (A.4) 



where rjB denotes one of the anomalous dimensions of the bosonic propagators under con- 
sideration, or ?7f in our case. The derivative dt acts on the k dependence of the cutoff 
function only (for an explicit representation of dt, see j3]). Some relations among the 
threshold functions are given by 



lf^if^''^{wF,WB,WB;r]i^,r]B) = /~!™^''^(wf,wb;'7^,'7b), 
/[™^''^(wf = 0, wb; v^, Vb) = 4™^''^(^f = 0, wb; v^p, Vb)- 



(A.5) 
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For our numerical computations, we use the linear cutoff functions proposed in [12] {y = 

r(y)= Q-l^^a-y), rF(y)= (^-^-l^e(l-y), (A.6) 

for which all integrals listed above can be performed analytically, yielding in the present 
context: 



~ (FB),d/ ^ N 

/™''^(wF,0;r/^,r/j.) 
J{FBB),d^^^^ wi, W2; V^, Tj^) 



1 



d — \ 1 + uip 
2 1 
2(1 + WyY 



2 



(A.7) 



1 _ + 



d + 2 ' (i + iy ' I+wfV^ d+1 



d{l+WFy{l+Wi){l+W2) 

1 



1 \ 



(i+2/ Vl+^f^F 



11 



d+1 



The representations of all other threshold functions for the linear cutoff can be looked up 
in Ha. 



B Nonperturbative running of the gauge coupling 

The infrared quantities serving as "physical observables" in the present work, such as the 
constituent quark mass or the eta boson mass, depend on the way we model the effective 
gauge coupling in the nonperturbative domain in our truncation. In order to gain more 
insight into this dependence, we study different gauge coupling f3 functions proposed in the 
literature in this appendix. Here we focus on theories within the QCD universality class 
which are sensitive to the infrared physics of the gauge sector. 

In Sect. EJ we used a (3 function with accurate two-loop behavior and an IR fixed-point 
at a* = 2.5. We denote this (3 function defined in Eq. serving as a reference as 
/?Ref in the following. Such (3 functions with a fixed point of the gauge coupling in the 
infrared have a long tradition in the literature and have frequently been discussed from a 
phenomenological viewpoint . Furthermore, some theoretical evidence for the existence 
of such a fixed point has been collected in certain nonperturbative approximation schemes. 
However, due to the lack of a unique nonperturbative definition of the gauge coupling and 
due to an inherent regularization scheme dependence of the /? function, a comparison of 
different theoretical approaches and a connection to phenomenology is difficult to make. 
Here we take a pragmatic point of view and use the various running couplings as effective 
ones which are implicitly defined by their use in our approach. 

Recently, an actual nonperturbative computation of the running coupling has been set 
up in the framework of truncated Schwinger-Dyson equations in Landau gauge [22] , reveal- 
ing an infrared fixed point; these results also receive some support from lattice calculations 
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[2n]- For our purposes, we use the representation given in for the running couphng, 



He + a,x^^+c,x^^y where «,SDE = 2.972, 



(B.8) 



and ai = 5.292GeV 



-2a2 



^2 



2.324, ci = 0.034GeV-2'=2^ ca = 3.169. This couphng is also 



normahzed to the standard value at the Z mass, and we identify x = A;^/(GeV)^. The /3 
function is given by Psde = ^tS'iDE- 

As a second example, we use the running coupling arising from a scheme called "Ana- 
lytic Perturbation Theory" j2Hl that has been devised for enforcing analyticity properties of 
the coupling in the time-like and space-like (Euclidean) region. For our numerical routine, 
we use the approximate (but two-loop accurate) representation 



fi'APT(^) 



(47r) 



1 



1 



kix) l-exp[/2(x)] 
In a; -|- ^ In x + An^, 



"^0 



(B.9) 
(B.IO) 



where we identify x = A;^/(1349MeV)^, so that this coupling is also normalized at the Z 



47r 
bo 



1.22 for 



mass. In the infrared /c — > 0, the coupling tends to the fixed point a=i.,APT 
A^c = 3 and A'f = 1. 

As a third example, we use a calculation of the running coupling based on a truncated 
flow equation that also revealed an infrared fixed point a*,FE jSH]- The corresponding 
/3fe function was obtained as an extensive multiple integral which we will not display 
here. Since this result holds for pure gauge theory, we incorporate one quark flavor in a 
"quenched" approximation by adding the fermionic part of the two-loop (3 function to the 
pure gauge result. This leads to an infrared fixed-point value of a^^pE — 3.43 it 0.01, where 
the theoretical error arises from an incompletely resolved color structure in jSH] • We would 
like to point out that the definition of the running coupling used in [21] agrees with the 
one of the present work. 

As a simple example for a running coupling which does not tend to an infrared fixed 
point, we employ a class of /5 functions that correspond to anomalous dimensions of the 
gauge field tjf which become constant for k ^ 0. This is realized by the choice 



-2 h 



9' 



167r2 



9^ 



(167r 



2^2 



exp 



'IGvr 



2\2 



fB.lll 



so that in fact rjF = Pri^./g'^ ?7* for A; ^ 0. For negative rj^,, the running coupling increases 
~ (l//c)l''*l for k 0. As explicit examples, we choose r]^ = —0.1 and t]^ = —0.5 for the 
numerical analysis. 

The results of the numerical integration of the flow equations are collected in Tab. El 
For the various (3 functions denoted in the flrst column, we listed the transition scale k^sB 
into the xSB regime and the generated fermion mass rrif in the next two columns. These 
results refer to calculations without axial anomaly and instanton-mediated interactions. 
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[3 function 


mf/MeV 


m^/MeV 




k^SB /MeV 


mf/MeV 


/?Rcf 


1765 


4438 


2.5 


423 


371 


/?SDE 


1777 


4226 


2.4 


457 


427 


/3apt 


563 


1990 


3.5 


4 


3 


/^FE 


903±2 


2117±1 


2.3 


243±2 


241±3 




513 


1793 


3.5 


11 


8 


/3j;.=-0.5 


1946 


4534 


2.3 


394 


395 



Table 1: Characteristic masses nif and for various nonperturbative (3 functions for the 
strong gauge couphng. The main uncertainty concerns the overall scale, whereas the ratio 
m^/rrif is relatively robust. We also show the scale of transition to x^B and the fermion 
mass rfif in absence of instanton effects. 

similarly to Sect.|3] In the last two columns the fermion mass with instanton contribution 
and the mass of the eta boson are given. 

Obviously, the quantities k^sB and mf in the calculation without axial anomaly are 
roughly correlated. Furthermore, Papt and /5r,,=_o.i lead to small values for k^sB and 
mf, since both approach larger values of the coupling only very slowly. The fermion and 
eta boson masses including the axial anomaly are not strictly correlated with the former 
quantities. The running of the gauge coupling enters these quantities over a wider range of 
scales, since sizable instanton contributions can already arise while the bound-state fixed 
point is still present. 

Nevertheless, our main observation is that the overall qualitative picture of the approach 
to xSB, even in the nonperturbative domain, is rather independent of the details of the 
gauge sector in our truncation. On the one hand, a strong gauge coupling g"^ > g"^ is all 
that is needed to trigger xSB; on the other hand, fermion decoupling cuts off any strong 
influence of the running coupling in the deep infrared. Quantitative results, of course, 
depend strongly on the flow of the coupling between the scale at which g'^ = g^ and the 
scale of fermion decoupling. This mainly concerns the overall scale, whereas mass ratios 
like rriri/mi turn out to be more robust. 

C Flow equations with axial anomaly in the broken 
regime 

Here we collect the flow equations for the various couplings in the broken regime, including 
the contributions arising from the axial anomaly. Let us begin with the scalar and fermion 
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anomalous dimensions: 

i2 _4 



V<l> 



(C.12) 



2C2(iVe)t;4 [(3 - mf^^^\Kh\ 0; r^^, r^p) - 3(1 - ml^''\KK\ 0; r/^, r/p) 

+t;4 /i^ [m[™^'^(K/i^ + 2kA<^; ?7^, t]^) + mS™^'^(K/;,^ 77^, r]^)] . (C.13) 

Including the appropriately adjusted fermion-boson translation as outlined in Sect. the 
flow equations for the minimum of the scalar potential and the scalar self-interaction read: 



~ (FB),4/ ,2 



dtK 
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4,(F),4^^^2., 

(/5C^?' + /5f>') (C.15) 
4Am / 



2r/^ A^ + 2v^ Aj r]^) + 94(^ + 2fi:A<^; r/^)] - 8iVcy4 4''^''(«:/i'; V^) 
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2v^ 
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The Yukawa coupling flows in the broken regime according to 

-8(3 + i)C2{N,)v, g^h^ /^?^''(/^/^^ 0; rj^, Vf), (C.16) 
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and the flow of the axial anomaly is given by 



2v^ jJV, 
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In these equations, the quantity is also modified 
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(C.17) 



(C.18) 
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whereas /3f remains the same. 

Act 
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